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Abstract. A nonequilibrium Green's functions (NEGF) approach for spatially 
inhomogeneous, strongly correlated artificial atoms is presented and applied to 
compute the time-dependent properties while starting from a (correlated) initial 
few-electron state at finite temperatures. In the regime of moderate to strong 
coupling, we consider the Kohn mode of a three-electron system in a parabolic 
confinement excited by a short pulsed classical laser field treated in dipole 
approximation. In particular, we numerically confirm that this mode is preserved 
within a conserving (e.g. Hartree-Fock or second Born) theory . 
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1. Introduction 

'Artificial atoms' (AA) are inhomogeneous quantum few-particle systems confined in a 
trapping potential and show bound (discrete) electronic states, as they are occurring in 
real atoms [T]. Most artificial atoms are realized in an (isotropic) parabolic confinement 
and quantum dots are a synonym convention for these systems, e.g. Refs [2J [3J. 
But AAs are also formed by ions in Penning and Paul traps, charge carriers in 
semiconductor heterostructures (quantum wells), or electrons in metal clusters. AAs 
in one- (ID) and two-dimensional (2D) entrapment show interesting properties far 
from ideal Fermi-gas behavior including ring and shell structures — see the 2D ground 
state configurations for N — 2, 3, . . . , 7 electrons displayed in Fig. [T] This is due to 
the externally controllable electron-electron interactions which, in particular, induce 
correlation phenomena. 

Previous theoretical investigations of AAs mainly concern ground state 
calculations and approaches to thermodynamic equilibrium, see Refs. [3J H] and 
references. The aim of this paper is to present a time-dependent theory including 
moderate and strong correlations. Following a quantum statistical approach, we 
thereby want to consider the system's response to a (strong) short-pulsed laser field 
while being embedded into an environment of finite temperature. To this end, in Sec. [3] 
and [I] we describe the NEGF method applied to compute the (correlated) initial state 
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(a) N = 2 




E/E* = 5.655 



(d) N = 5 




EjEl = 33.474 



(b) N = 3 




E/E* = 12.571 




E/E* = 46.792 



(c) N = 4 




E/E£ = 21.913 



(f ) N = 7 




EjE' a = 61.961 



Figure 1. Examples for the one-electron ground state (GS) density n^j s (r) of 
TV = 2, 3, . . . , 7 electrons in an isotropic 2D trap for /3 = 100 and A = 5.0. 
Figures (a) to (f) display the result of (symmetry broken) unrestricted Hartree- 
Fock calculations. The side length of the density plots measures 8x^. 



which is then propagated in time according to the Kcldysh/Kadanoff-Bayni equations. 
In Sec. [5] we present results for a three-electron A A in a ID and 2D trap geometry. 

2. Model 



The (i-dimensional A-electron Hamiltonian of the artificial atom in a time-dependent 
laser field E(i) described in dipole approximation reads 



m) = E - 



2m! 



^ + eE(i)r 4 + £ 



N 



i<3 



Aire Ti 



(1) 



where the effective electron mass is given by to*, the frequency luq adjusts the 
confinement strength (and hence the density in the AA), e is the elementary charge 
and e denotes the background dielectric constant. Further, the d-dimcnsional electron 
coordinates originate from the trap centre and ry = |r< — rj\. 

Hamiltonian can be rewritten in dimensionless form, using the transformations 
{E — * E/Eq,t — > t/xq}, where Eq — huio is the confinement energy and the 
(oscillator) length Xq = h/ {m* e ujo) denotes the characteristic spatial one-electron 
extension in the AA. Consequently, the system is characterized by a single coupling 
(or Wigner) parameter A, which relates the characteristics Coulomb energy Ec — 
e 2 /{A-kexq) to Eq according to 
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with as being the effective electron Bohr radius. The Hamiltonian ([!]) then transforms 
into dimensionless form 

N N 

HxW = g E (- V ' + r i + + X E — ' ( 3 ) 

i=l t<j lJ 

with 7(t) = eE{t) / yJhm* e ujQ. The coupling parameter A adjusts the influence of 
electron-electron interactions and (quantum) correlations: In the case of A <C 1, the 
artificial atom behaves similar to an ideal Fermi gas. For A ~ 1, the equilibrium 
state of the AA is Fermi liquid-like, whereas in the limit A — > oo, it is Xq 
as, and quantum effects vanish in favor of classical interaction dominated charge 
carriers [10]. For moderate coupling (A > 1) the AAs typically show spatially well 
localized carrier density including Wigner molecule (Wigner crystal- like) behavior [3], 
see Fig. [I] Moreover, if the AA is not in its ground state (GS), one has to take into 
account thermodynamic fluctuations due to a surrounding heat bath of dimensionless 
temperature /3 _1 = fc^T '/-Eg. Below, all presented results arc related to the system of 
units {xq , Eq } and in the definition of the NEGFs we take H = 1 . 



3. Preparation of equilibrium states 

Introducing electron annihilation (creation) operators ip^\riti) acting in the 

Heisenberg picture at a space-time point 1 = (ri,£i), the second-quantized form of 
d3l reads 



H x (t 1 )= I d<Wt(l)fr°(lMl) (4) 

d d rd d f ft(i)ft(T) w(n - ri) 4>(1) i>(i) , 



with the one-electron energy h°(l) = (— + rf)/2 + 7(fi)ri and the interaction 
w(ri — rj) = A |ri — rj| _1 . In the following, we study Hamiltonian Q at finite 
temperatures by means of the one-particle nonequilibrium Green's function 

G(l, 1), which is defined on the Schwinger/Keldysh contour C (see e.g. Refs. [SUHj) as 

G(M) = -i<T c V>(l)^t(i)} j (5) 

where Tq denotes time-ordering on C. G(l, 1) obeys the two-time Keldysh/Kadanoff- 
Baym equation (KBE)[7] 

[»a tl -ft(l)]G(l,l) = <Jc(l-l)-»j^2W(l-2)G 12 (l,2;l,2+) , (6) 

and its adjoint, where W(l — 1) = 8c{t\ — ti)w{r\ — rj) acts instantaneously with 
contour delta function 8c, and 2 + indicates the time limit £2 ~^ £2 + 0. Here, the 
two-particle Green's function, 

G 12 (l,2;l,2) = (-if(T e ^t(i)^t(i)^t (2) ft (2)) , (7) 

appears as a generalization of the two-particle density matrix. In terms of 
G12, we can formulate all relevant many-body approximations: For instance, 
substituting Gi 2 (l,2;l,2) -► G(1,1)G(2,2) - G(1,2)G(2,1) yields the Hartree-Fock 
(HF) approximation. Second and higher order approximations (e.g. second Born or 
GW) , can be systematically obtained by diagram expansions known from ground state 
many-body theory and are valid for both equilibrium and nonequilibrium situations. 
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When in Eq. Q, we consider 7(ii) = for — oo < t\ < to, the AA stays in 
thermodynamic equilibrium until at a time t\ > t j(ti) becomes nonzero. Without 
loss of generality we thereby can take to = 0. Specifying the time-independent single- 
electron part of Eq. (jij) as h°(ri) = /i(l)| 7 ( tl ) =0 , the KBE then reduce, for t 1 j < 0, 
to the Dyson equation 

[-^-/^(n^W^r) (8) 

= S(r)+Jd 2 fJ^ dfSf (n,f;r-f) G M (f,r i; f) . 

Here, the Matsubara Green's function being defined as G M (1,1) = G (ri,rj;r) = 
G(ri0 — iri,rj0 — ztj), with r = t\ — tj e [—/?,/?], characterizes the equilibrium 
(initial) state of the AA. Further, on the right hand side, we have introduced the 
self-energy Sf (1,1) = Sf (n,r i; r) according to -iW(l - 2)G 12 (1, 2; 1, 2+)\ h ^ = 
£[G](1, 2)G(2, 1). A conserving many-body approximation [7 , i.e. an approximation 
for £ that preserves density (continuity equation), total energy and momentum, can 
now be formulated in terms of a functional $ such that S(l, I) = (5<I>[G]/<5G(I, 1). 

At a given temperature /3" 1 , most of the equilibrium properties of the AA system 
(IT]), e.g. total energy, one-par ticle density and energy spectrum, are contained in 

o 4- Vi i2i -PV~iT»T-m line m Qj^ir» /I O o n r\ 4- o l^d 4- Vi d 1 itn 1 4- + j i. O T7V\r- 4- Vi a niimfln r>o 1 



G , see the formulas in Sec. 4.2 and take the limit t\ — > 0. For the numerical 
techniques applicable in solving Eq. ^ in matrix form see e.g. Refs. j3j[8j[9]. In HF 
approximation, the self-consistent solution can be written as 

n b — l 

G M (ri,ri;r) = j ^(n) ^(ri)fll(r) , (9) 
m=0 

= //3(e m - m) e^"^ - e-^-")/^"") + 1) , 

with interaction renormalized (effective single-electron) HF orbitals m (r)[lO , 
quantum numbers m = 0, . . . , n& — 1, discrete energies e m , and a chemical potential /i. 
Beyond HF level, G M will be no longer diagonal in the functions <j) m , and the respective 
occupation probabilities will deviate from a Fermi-Dirac distribution fp{t m ~ t 1 ) due 
to additional electron scattering processes. In particular, the inclusion of electron- 
electron correlations leads to orbital-dependent energy shifts and broadening in the 
HF spectrum a(cu) = J2 m — e ™), see Ref. [3]. 



4. Time-propagation of initial states 

When for t\ > the laser field is switched on and 7(^1) ^ 0, the quantum state of the 
AA evolves in real time according to the KBE, Eq. ^ and its adjoint. Thereby, being 
computed from the Dyson equation (Jsj) in a self-consistent manner, the Matsubara 
Green's function serves as initial (Kubo-Martin-Schwinger) condition for the time- 
propagation. In particular, for to = 0, one has 

GfaO-iTu^O-iTs) =i[G M (r 1 ,r i ;ri)-G M (r 1 ,r I ;-Ti)] . (10) 

Beyond mean field level, all relevant initial correlations are taken into account via G M 
and, consequently, evolve in time, leading to a correlated TV-particle dynamics. 
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4-1. Solving the Keldysh-Kadanoff Baym equations 

The expansion of the NEGF in terms of a HF basis, see Eq. (|9j, advises us also to 
solve the real-time KBE in matrix formJSlQT]. This means that we generally consider 

"6 — 1 

G(l,l) = 4>* m {*l)<t>n{n)9rnn{tl,tl) , (H) 

m,n— 

with time arguments t±,ti on the contour C, the coefficient matrix g m n{ti,ti) = 
S(ti, ii)5 l m„(*i, *i) — 0{Hiti)9mn{^iiH) °f dimension n b x n^, and steady-state HF 
orbitals </> m (r) which generate a complete orthonormal set. Hence, g m nitiiti) = 
—i(Tc c m (ti) cjj(tj)) are just the NEGFs with respect to the operators & m (c m ) that 
create (annihilate) an electron in the state m. Consequently, the diagonal elements 
[to = n] are directly related to the occupation numbers (probabilities) of the HF 



orbitals, cf. Eq. (IT), whereas the off-diagonal elements [to ^ n] are connected with 
'interband' excitations, i.e. the transition probabilities between the energy levels. 

Inserting expression (111, Eq. ([6| and its adjoint then transform into integro- 
differential equations for the matrix elements g m n(ti,ti) = (g)mn(iijii) : 

[id tl -h(ti)]g(ti,t I ) =5 c (ti-ti) + Jd£H{ti,i)e(I,ti) , (12) 

g(*i,tO [-id tl -h(ti)] =*c(*i-ti) + j^dtg(ti, , (13) 

where the single-carrier energy is given by h mn (ti) = J d d ri 0^ l (r 1 )/i(l)0„(r 1 ), and 
in the adjoint equation ( |T3| ) the operators are acting to the left. The explicit form of 
the self-energy matrix E nm (*i.*i) at the HF level is 

n b — l 

S mn(*l)*l)= -i^Scih-ti) ^ [Wmn,kl~ W mlM]9kl{ t lM) > ( i4 ) 

k,l=0 

with the two-electron integrals w mn ,ki defined as 

d d rd d f ^(r) ^(f ) u;(r - f) 0„(r) 0,(r) . (15) 



The detailed structure of E(ti,*i) = S HF (ti, t T ) + S corr (ti, t{) in second (order) Born 
approximation is given e.g. in Refs. [3j ITT| IT2j. 

For the real-time arguments in g(ti, tj), it is useful to introduce relative and centre 
of mass (cm.) variables, t c , m , = (t\ + £j) /2 and t Te \. =t± — t\. The Green's functions 
with respect to the cm. time then account for the statistical (thermodynamic) 



properties of the artificial atom, see the definitions (16 1 to (19 1, while quantities with 



respect to t Te \. carry the dynamical (spectral) information, cf. a{uS) as defined at the 



end of Sec. 4.2 Moreover, we note that Eqs. ( pL2[ ) and (13 1 can, in HF approximation, 
be further simplified to a single-time equation involving only t c m . More technical 
and numerical details for the time-propagation of the Green's function matrix are to 
be found in Refs. 0QI]. 

4-2. Dynamical quantities 

The spatial one-particle density in the AA and the HF orbital-resolved occupation 
probability of state to are 

"6 — 1 

(n)(r 1 ,i 1 )= -iG{l,l + ) = -i O-i)0n(ri)Smn(*i.*i)> (W) 

rn.n— 
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HF orbital energy e m [wo] 

Figure 2. Thermodynamic properties of TV = 3 charge carriers in the ID AA 
with A = 2.0 and (5 = 2. a) HF distribution function /a(e m — A 1 ) (solid curve) 
and renormalization at the second Born level (dashed curve). The triangle marks 
the position of the chemical potential in HF approximation. Further, the ideal 
energies = TO + 1/2 (marked with crosses) for A = are shifted by the HF 
value of fj,. b) Density profile (n)(x) for the ideal system as well as for the HF 
and second Born approximation. 



(n m )(ti) = g< m (tx,h) . (17) 

In a conserving approximation, the electron number N(t±) = ^2 m (n m ) (t\) in the AA 
is preserved in agreement with the continuity equation, i.e. 9 tl (n)(l) +div(j)(l) = 
with current density (j)(l) = -|{[V rj - V rj ] G < (r 1 t 1 , r I i I )} ri=ri . 

The relevant energies involved are the kinetic and potential energy 

n b -l 

(-Bkin+pot)(*l) = ~i £ [*mn + «mn(*l)]5nm(*l)*l) ! ( 18 ) 
m,n— 

with corresponding definitions t mn = J d d ri (/^(ri) [— V, /2] 4> n { v i) and u mn (ii) = 

Jd d n 4r m (n) [rf/2 + 7(*i)ri]<k(n). thc HF ener sy 

. "6-1 

(£ H F>(ti) = - \ £ S^(fi) 9 < m (ti,ti), (19) 

m,n— 

and the correlation energy (E COI r)(ti) = -| / c di Tr{S corr (ti, t) g(i, i+)}. 

Other interesting quantities are e.g. the time-dependent dipolc moment 
(d)(ti) = —ief d 2 ri{vi G < (viti 1 rjii)} ri=ri or the one-particle spectral function[3] 
a(u,t c . m .) = iT, m I dt rei. e iut " l \g> m (ti,ti) -5< ro (ti,*i)]. 

5. Numerical results 

In this section, we study the dynamical properties of a three-electron AA, when, 
initially in thermodynamic equilibrium, the system is excited by a single few-cycle 
laser pulse described in dipole approximation, cf. Eq. More precisely, the field is 
linearly polarized in cc-direction and has the time-dependence 

7 (t) = So e -(*-' ! ) 2 /(2rf) cos(w;(i _ t[)) 5 (20) 



Nonequilibrium properties of strongly correlated artificial atoms 



7 




10 20 30 40 50 60 
time t [^q" 1 ] 



10 20 30 40 50 60 70 
time t [^o" 1 ] 



10 20 30 40 50 60 
time t [uJq 1 ] 



Figure 3. Nonequilibrium behavior of the TV = 3 AA system (ID) with 
parameters A = 2.0 and /3 = 2. a) to c) show the mean-field dynamics of the 
HF orbital occupation probabilities (n;) (t) at a near-resonant laser frequency 
ldi = 1.2hujQ and three different pulse durations tj. The corresponding pulse 
shapes -y(t) are indicated above the figures. 



where £q — £i/^/2tt denotes the amplitude of the electric field, the Gaussian envelope is 
centred at ti , the pulse duration (variance) is given by r/ , and the oscillation frequency 
is uji, cf. Figs. |3]a)-c). 

As the response characteristics of the quantum system we define the amount 

of energy that has been absorbed from the laser field for a fixed frequency oji, i.e. 

S(wi) = {EtatUtt ->«>)- (£tot)(0) . (21) 

This quantity together with the time-dependent occupation probabilities (n m )(t) 
allows for the determination of (off)resonant nonequilibrium behavior (and nonlinear 



effects), see Sees. 5.1 and 5.2 

The nonequilibrium behavior of the quantum system is theoretically well 
known: Driven by the laser field, the AA exactly responds according to the excitation 
of the centre of mass (Kohn or sloshing) mode. This is obtained from the Kohn 
theorem, and its generalization to the case of an additional external dipole field, see 
Ref. [13 and references therein. Its statement is that, independent of dimensionality, 
the centre of mass coordinate R(t) = N^ 1 J2i r iW °f a parabolically confined, 
interacting electron system performs (equivalently to a single particle in the AA) 
the motion of a forced harmonic oscillator, |R| + Wq|R| = NeE(t)/m*. Furthermore, 
this effect is accompanied by a rigid translation of the density profile (n) (r) , since the 
particle interaction appears only in the relative Hamiltonian and [-ff c .m. , -ffrci.] = 0. 
The key point in the present study is, however, that the Kohn theorem also holds 
when the interaction is treated approximately, as long as density, total energy and 
momentum are preserved (conserving approximation), and also applies to zero and 
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<«Hp)(t) 



(Epo.}(«) 



(JW(!) 



(-Eto.Xf) 
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10 20 30 40 50 60 10 20 30 40 50 60 70 10 20 30 40 50 



time t u)j 



ii 



time i \uj, 



-ii 



time i U, 



-ii 



Figure 4. For the laser irradiated artificial atom [ojj = 1.25uo] also considered in 
Fig. [3] a) to c) show the mean-field dynamics of the relevant energy contributions 
in dependence of the pulse duration. For tj = 10, the AA returns after the 
excitation close to its initial state (off-resonant situation). Whereas, in Figs, b) 
and c) the spectral width of the laser frequency is essentially increased compared 
to a), leading to the AA remaining in an excited state of the Kohn mode (resonant 
case). Potential and kinetic energy thereby oscillate out of phase with exactly 
double confinement frequency luq, while the HF energy stays constant. 



finite temperatures — for the proof see Ref. [T3] ■ 

The following mean-field results for ID and 2D have been obtained from NEGF 
calculations with up to rib = 40 HF orbitals. The main limitations of the approach 
are thereby (i) the basis size, which sets the dimension of the time-evolution matrix 
U(£i) = exp(— i [h(ti) + S HF (£i)] ti) to be computed in each time-step (diagonalizing 



h + S HF ), and (ii) the two-electron integrals Wij^i, Eq. (151, that generally require 
large memory resources [scaling with nf] and need to be processed very frequently 
in the self-energy expression S(ti,tj). With more than 50000 time-steps needed 
to achieve convergence, this results in computing times of typically several hours 
on a single machine. For the correlated time-evolution of the AA in second Born 
approximation, the propagation must be carried out in the whole two-time plane 
(ti,tj). This is an even more intricate task as one needs to compute all the higher 
order collision integrals on the r.h.s. of the KBE. However, these calculations are 
currently near completion — examples are to be found in Ref. [8] and for applications 
on real atoms and small molecules see Refs. [51 ITT1 |T2~] . 

5.1. ID case 



For the ID AA calculations in equilibrium and nonequilibrium, we, in Eq. (15 1 have 
replaced the pure Coulomb interaction w(x— x) by A [(x — x) 2 +a 2 ]^ 1 / 2 with a = 0.1 as 
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Figure 5. Response characteristics 3(u>i), Eq. \21\ , for the pulsed laser excitation 
of N = 3 electrons in the ID artificial atom. The system parameters are as 
in Fig. [2] The AA shows resonance behavior at the confinement frequency 
only, i.e. for ui\ = u)q, and responds via the cm. motion (Kohn mode) — rigid 
translation of the whole density. With increasing pulse durations t; (sharpened 
laser frequency) the resonance curves become more and more peaked. 



a regularization parameter. This is necessary to make the integrals w mn M finite and, 
in a physical interpretation, allows for a small transversal spread of the one-electron 
wave functions [2, . In 2D, we used a = as the integrals converge. 

We consider three electrons in a ID artificial atom at (3 = 2. With A = 2.0 
the system is tuned into the crossover regime between Fermi liquid-like and crystal- 
like behavior. The equilibrium properties can be read from Fig. g] for the HF and 
second Born approximation. For the HF energy spectrum and the corresponding 
distribution function fp(e rn — /i) of the equilibrium (initial) state, including its 
collisional renormalization in second Born approximation, see Fig. |] a). The one- 
electron density (h)(x) is displayed in Fig. [2]b). Compared to the HF result (solid 
curve, (£^tot) = 8.790), here, the inclusion of electron-electron scattering (dashed 
curve, (i? t 2 " t dB ) = 8.941), leads to a considerable reduction of the density modulation, 
which is accompanied by an increase of the total energy of about 1.7 %. 

Starting from the HF Green's function G M (ri, rj; r), the AA was now propagated 
in time under the presence of a laser field (centred at i; = 25) with amplitude £1 = 0.5 
and frequency u>i = 1.25wo- What happens to the orbital occupations and the energies 
for different pulse durations 77 is shown in Figs. [3] and [3] In all cases, gradually, the 
HF orbitals m < 3 become depopulated and the states m > 3 analogously become 
populated with preservation of N(t). Oscillations of the increased total, kinetic and 
potential energy and (n m )(t) thereby occur with twice the confinement frequency. 
In Figs. [3] a) and [4] a), respectively, the laser excitation is such that the iV-particle 
dynamics is decelerated and almost freezed after the pulse has passed. Consequently, 
we nearly recover the initial state characterized by (n m )(0) and (iS to t)(0)- Also, for 
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Figure 6. Thermodynamic initial state of the 2D artificial atom (N = 3) at 
A = 2.0 and /3 = 2 (in HF approximation), a) Single-electron density profile 
(n)(r) which is rotationally symmetric. Figs, b) to g): Energetically lowest 
spatial, unrestricted HF states tp m (r) with orbital energies eo>---; e 5 where e\ 
and £2 as well as eg and £4 are degenerate). The arro ws in c) and d) mark the 
direction of polarization of the laser field 7(t), Eq. 1201 ■ 




occupation probability (h m ) (t) 



1 pulse 



Figure 7. Mean-field dynamics of the HF orbital occupation probabilities 
(ftm) (t) = n m for three charge carriers in a 2D artificial atom with coupling 
parameter A = 2.0 and inverse temperature /3 = 2. The laser frequency is again 
near-resonant, lui = f.25<^o, and the pulse duration is t; = 5. In the initial state 
the occupation numbers ni and ri2, n$ and as well as ng and n-j are practically 
pairwise degenerate, compare with the energy spectrum displayed in Fig. [6] 



different pulse durations, i.e. different spectral profiles of the laser, the maximum 
laser energy absorption is observed at the confinement frequency loq, cf. the response 
function S(<x>i) in Fig. [5] In addition, the single resonance-peak at ecu — luq sharpens 
with the increase of 77. Consider now the spatial dynamics. We observe in all cases, 
that the center of mass of the AA R(i) performs a harmonic oscillation with frequency 
too while the whole density profile itself is translated rigidly. Accompanying this fact, 
.Ehf is constant in time, see Fig. [4] Thus, we numerically confirm that the Kohn 
theorem is satisfied. 
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5.2. 2D case 

For the three-electron AA in 2D, we have chosen the same system parameters, 
A = 2.0 and /3 = 2 — however, no regularization parameter a was needed in the 
two-electron integrals w mn ,ki- The (almost) rotationally symmetric density profile 
for the equilibrium state is shown in Fig. [6] and indicates a ring-like structure. It is 
instructive to note, that the unrestricted HF solution of the Dyson equation leads to 
orbitals 4> m (r), that are in general arbitrarily oriented in space, compare Fig. [6] b) 
to g). Together with the energetically degenerate states m = 1 and 2 (to = 3 and 
4, etc.), this has the following consequence on the dipole excitation: As degenerate 
orbitals can be differently oriented relative to the laser field, in the time-evolution of 
the artificial atom this degeneracy is lifted. In the present case, the orbitals to = 1 
and to — 2 (see Fig.[6]c) and d)) arc almost aligned with the diagonals (dotted lines), 
nevertheless the small deviations are sufficient to clearly influence the evolution of the 
occupation numbers (n m )(t), see (hi) and (712) as well as (713) and (n.4) in Fig. [7] 

6. Conclusion and outlook 

We have presented an analysis of femtosecond relaxation of few-particle quantum 
dots during and after a short laser pulse. The method of NEGF wa shown to be 
efficient to describe the dynamics even in the range of strong Coulomb correlations. 
Numerically, the Cm. mode excitation can serve as a very sensitive test for the 
NEGF calculation [and any other numerical code] involving quantum many-body 
approximations Ref. [T3] , 
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